Origin of magnetism and quasiparticles properties in Cr-doped Ti02 
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Combining LSDA+C/ and an analysis of superexchange interactions beyond DFT, we describe the 
magnetic ground states in rutile and anatase Cr-doped Ti02. In parallel, we correct our LSDA+C/ 
ground state through GW corrections (GVF@LSDA-|-C/) that reproduce the position of impurity 
states and the band gaps in satisfying agreement with experiments. Because of the different topo- 
logical coordinations of Cr-Cr bonds in the ground states of rutile and anatase, superexchange 
interactions induce either ferromagnetic or antiferromagnetic couplings of Cr ions. In Cr-doped 
anatase, this interaction leads to a new mechanism which stabilizes a ferromagnetic ground state, 
in keeping with experimental evidence, without the need to invoke F-center exchange. 

PACS numbers: 71.10.-w,75.47.Lx,71.15.Mb,75.30.Et 



Understanding the physics of pure and doped 
transition- metal oxides (TMO), in order to have reliable 
predictions of their ground and excited states, represents 
a major challenge in fundamental and applied research 
today. The main difficulty for an ab-initio description 
of the TMO wavefunctions stems from the partly local- 
ized character of their 3c?-electron states, which often 
leads to deep interplays of charge, spin and orbital de- 
grees of freedom [l[ . For this reason, Density Functional 
Theory (DFT), within its traditional functionals, fails to 
describe their electronic and magnetic properties 0, as 
Kohn-Sham wave functions are too delocalized. In order 
to cure such a drawback, alternative schemes like hy- 
brid functionals [3, SIC [j] and LDA-hJ7 % have been 
used in recent years. The latter, in particular, due to 
the strong on-site Coulomb repulsion, leads to a more 
pronounced localization of 3c? wave- functions that bet- 
ter represents their behaviour and can therefore be used 
as a zeroth-order for subsequent many-body corrections 
{GW@LST)k+U approach (l-Q). 

The recent discovery of room temperature ferromag- 
netism (FM) in TM-doped Ti02 ^ has triggered an enor- 
mous experimental and theoretical effort to understand 
the magnetic properties of a class of compounds (Ti02, 
ZnO, SnO) that show correlation effects when doped with 
TM impurities. Though F-center exchange is often in- 
voked to justify FM in these and analogous systems Q, 
recent ab-initio calculations l3, III point to its exclusion, 
so that the actual physical mechanism remains an open 
question |12l4l7| . Theoretical results with advanced func- 



tionals on TM-doped Ti02 [18|, [19| allow to qualitatively 



reproduce band gaps, but do not explain FM. 

In this work we describe magnetic ground states quan- 
titatively by evaluating superexchange (SE) corrections 
to the LSDA-I-J7 zeroth-order magnetic energies and dis- 



cuss their relevance for the stability of magnetic ground 
states in connection with the topology of the host lat- 
tice. Our approach goes beyond DFT (that does not 
consider dynamic fluctuations) and current implementa- 
tions of Dynamical Mean-Field Theory, limited to atomic 
dynamical correlations [20]. The correctness of the 
zeroth-order LSDA-I-C/ calculations is confirmed by our 
GW^@LSDA-|-?7 quasiparticle results, reproducing pho- 
toemission experiments with satisfying agreement. With 
this calculation scheme, we perform a systematic analy- 
sis of rutile and anatase Cr-doped Ti02 and confirm that 
F-center exchange is not active in Cr-doped Ti02, as O 
vacancies are F+^-centers. Moreover, the investigation of 
SE interactions highlights a strong interplay between the 
topology of the local structure, the presence of an oxy- 
gen vacancy and the orbital occupation of impurity lev- 
els, leading to the conclusion that such oxygen vacancies 
can boost FM via a different physical mechanism than 
F-center exchange. The influence of oxygen vacancies on 
the positions of impurity states is also investigated. Mu- 
tatis mutandis, our approach can be adapted to a large 
class of doped TMO and it paves the way for a critical re- 
examination of magnetism in all those compounds where 
t2g levels play a major role and a simplistic application 
of the Goodenough rules [2l| is not appropriate. 

Calculation framework. Calculations are performed 
within the PAW scheme using the ABINIT package [22 1. 
Ti and Cr semicore states are included in the valence. 
We used supercells with 48 atoms sampling the Brillouin 
zone with 2x2x3 and 2x2x2 Monkhorst-Pack k-mesh 
for rutile and anatase. Structural degrees of freedom are 
optimized by minimizing quantum mechanical forces and 
total energy. GqWq calculations are perforrned using the 
Godby-Needs plasmon-pole approximation ^23] , and with 
2000 bands for the Green's function. SE calculations are 
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TABLE L Band gaps of majority and minority spins for ru- 
tile (r) and anatase (a) in LSDA+[/ and GoWo@'LSDA+U . 
Positions of impurity states of given spin polarization with re- 
spect to the top of the valence band are shown in parenthesis. 
In bold, values to be compared with experimental data (see 
text). All values are in eV. 



LSDA+C/ (t; i) 



GoWo it, ;) 



ICr (r) 
2Cr (r) 
2Cr-0„ 
ICr (a) 
2Cr (a) 
2Cr-0„ 



(r) 



2.0 (0.6); 2.0 
1.8 (0.5); 1.7 

2.2 (1.5); 2.2 

2.3 (1.0); 2.4 
0.5; 1.8 
1.6; 1.6 



1.9 (0.9); 1.9 (1.1) 
1.3 (0.6); 1.0 

2.0 (0.2); 1.9 
2.7 (1.1); 2.7 (1.4) 
2.4 (1.2); 2.5 (1.3) 

2.5 (0.7); 2.7 



described below. 

At first, we confirm that LSDA is a poor zeroth-order 
approach for TMO: a direct application of GqWq on the 
LSDA half-metallic state does not open a gap and the 
only effect of the GqWq self energy is to push down va- 
lence states with respect to the conduction band [see Fig. 
[Ija)]. Since the states near the Fermi level are largely of 
t2g character, we describe their partly localized behaviour 
by hSBA+U. We consider U=5.0 (5.2) eV for Ti (Cr) 



3d states, from the literature [2J], and use the values of 
Hund's couphng J as obtained from LSDA: 0.80 (0.84) 
eV for Ti (Cr), as J is almost unsensitive to screening 
effects. The LSDA-hC/ and GoWo@LSDA-)-C/ density of 
states (DOS) for the cases of one Cr substitution, two 
Cr substitutions, and two Cr in combination with an O 
vacancy are shown in Fig. [T]for the most stable configu- 
rations both in rutile and anatase. We also show the FM 
case for anatase with 2 Cr and 1 O vacancy. Our cal- 
culations correspond to Cr concentrations of 6.25% and 
12.5%. The band gaps and positions of impurity states 
are summarized in Table L 

LSDA+U results. The LSDA+[/ functional opens a 
gap at the Fermi energy for both rutile and anatase 
phases (with and without O vacancy). The band gaps 
and positions of impurity states are summarized in Table 
I and will be compared with available experiments later 
in the text. Here we focus ourselves on discussing the in- 
terplay between local structural, electronic and magnetic 
properties. The DOS of the impurity states is mainly 
composed of Cr t2g states. Cr-doping favors the forma- 
tion of an O vacancy, as the formation energy is lowered 
with respect to the undoped cases by 1.2 and 1.6 eV for 
rutile and anatase, respectively. For the lattices with 2 
Cr ions, such an O vacancy is found to be more stable at 
a bridging position between the two impurities, in both 
lattices. Interestingly, the most stable configurations in 
rutile and anatase, both with and without oxygen va- 
cancies, are characterized by a topological difference in 
the relative position of the two Cr ions: nearest neigh- 
bors (NN) in rutile, connected along the tetragonal axis 
(bond-angle Cr-O-Cr ~ 90°) and next nearest neighbors 



(NNN) in anatase, connected along a direction perpen- 
dicular to the tetragonal axis (bond angle ~ 180°). This 
result is also valid for LSDA calculations and is probably 
favoured by the relative closeness of NNN in the anatase 
lattice (3.8 A) compared to rutile (4.6 A). Furthermore, 
we found that the electron charge density at the O va- 
cancy is negligible (see Supplementary Material), i.e., the 
vacancy is a F++-center 25], thereby suggesting that F- 
center exchange [9] is not active in this case and that the 
origin of magnetism should be looked for elsewhere. 

Spin polarization mainly occurs at the Cr site (~ 99%) . 
Without oxygen vacancies, we obtain a total magnetic 
moment of 2 /is/supercell, i.e., the two t2g electrons in a 
triplet configuration. Removal of one neutral O atom in 
the supercell leaves two extra electrons that reduce two 
Cr'^^ to Cr^"*" ions leading to a bigger average magnetic 
moment of 3 /is/Cr atom, in keeping with the values 
2.6-2.9 hb/Ct atom observed in Refs. 13:|26[. We also 



performed fully relativistic non-collinear LSDA-|-t/-|-spin 
orbit (SO) calculations to check whether orbital angular 
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momentum, L, could have been unquenched by SO inter- 
action in t2g orbitals, but we found no relevant changes 
in the magnetic moment (L ~ 0). 

In Fig. m we plot the total LSDA+C/ energy for the 
FM and the AFM configurations as a function of the in- 
terdistance of the two Cr impurities. All most stable con- 
figurations are FM, except for anatase with O vacancy, 
that is AFM (but nearly degenerate with the FM case: 
AE^12 meV). 



3.0 3.5 4.0 4.5 5.0 5.5 6.0 6.5 




3.0 3.5 4.0 4.5 5.0 5.5 6.0 6.5 
Cr-Ct distance [A] 



3,0 3,3 3.6 3.9 4,2 4,5 4,8 5,1 
Cr-Cr distance [A] 



FIG. 2: (Color online) FM (filled circles) and AFM (open 
squares) LSDA+L'' total energy vs interdistance of the two 
Cr impurities. The energy of the most stable configuration is 
set to zero. All values are in eV. 

SE contributions to magnetic ground state. We an- 
alyze here the stabilization of the FM ordering due to 
non-local spin-orbital correlations that are not included 
in the LSDA+J7 starting point. It has to be noted that 
exchange interactions for Cr t2g electrons do not neces- 
sarily follow AFM Goodenough SE rules t21j , like in the 
case of Mn-eg electrons in perovskites. The sign of ex- 
change interactions is affected by the character of Cr 
electrons (contrary to Cg electrons, they do not point to- 
wards ligand anions) and by the different connection of 
Cr ions in the two lattices, i.e., a ^ 90° TM ion-anion-TM 
ion configuration in rutile, and a ^ 180° TM ion-anion- 
TM ion configuration in anatase. In what follows, we 
provide the sign and we estimate the expected order of 
magnitude of the leading terms of SE interaction ener- 
gies for the four configurations corresponding to Fig. [21 
We find that only rutile with O vacancy has an AFM ex- 
change energy, whose magnitude is however insufficient 
to overcome the FM LSDA+U stabilization energy. All 
other three cases are instead characterized by FM correc- 
tions. In particular, the FM SE energy of anatase with O 



vacancy overcomes the 12 meV of the AFM LSDA+[/ so- 
lution and leads to a FM ground state, as experimentally 
found. 



Consider rutile at first: we refer to the model of |27l.l28| 
for exchange of t2g electrons with a 90° TM ion-anion- 
TM ion bond-angle. The exchange determined by the 
direct overlap of in-plane t2g orbitals (e.g., dzx orbitals 
in the zx plane) is the leading term both with and with- 
out O vacancies. In the former case, orbital degrees of 
freedom of Cr'^+ ion are frozen and the magnetic interac- 
tion is AFM, at odds with the case of 90° TM ion-anion- 
TM ion of Cg electrons, which is FM 29]. The magni- 
tude of this AFM exchange energy can be estimated as 
^AFM ^ t^/(U2 + J) -22 meV, where t ~ 0.36 eV is 
the hopping amplitude of the 2 Cr ions [l^l and U2 — U 
[3^ . This AFM correction is however not sufficient to 
overcome the FM LSDA+[/ contribution shown in Fig. 
2(c). The case without O vacancies is different because 
Cr"*"*" ions are characterized by one empty orbital and, 
because of the Coulomb repulsion, 2 Cr ions facing in 
the xy plane have just one dzx orbital filled. This config- 
uration favours FM-coupling , with an energy gain of: 
£:^*^ ~ {J/U2)t^/{U2 - J) ~ 5 meV, slightly increasing 
the FM energy gap of Fig. 2(a). 

For anatase, the change in topology (~ 180° TM ion- 
anion-TM ion configuration) also determines a change in 
the physical mechanism leading to exchange interactions. 
The most interesting situation is that of anatase with 
O vacancies, as AFM and FM LSDA-I-J7 solutions are 
nearly degenerate. Consider the hopping process around 
the O vacancy, through the t2g-empty NN Ti-site. This 
is qualitatively described in Fig. [3] an electron from 
Cri hops to Ti2; then from Ctt, to Cri; then from Ti2 
to Crs; the opposite path is also possible. This process 
is allowed only in case of FM coupling of Cr'^^ ions and 
is in competition with the two AFM processes where a) 
both electrons hop to Ti2 with opposite spin; b) the two 
electrons directly hop from one Cr'^^ ion to the other, 
as for rutile, but with a much smaller hopping ampli- 
tude, ii3, as ions are much further. Notice that second- 
order hopping from Cr ions to NN Ti ion does not lead 
to effective exchange, as the energy gain does not de- 
pend on the magnetic coupling. Quantitatively all these 
processes can be described through a Hamiltonian anal- 
ogous to that of Ref. (s^l, with intra-orbital Coulomb 
repulsion t/^' at Ti-sites. We get a FM SE energy of: 



E 



SE 
FM 



= ^tH 
E' 



-^AFM 



AFM 



is/Act 

4, 



"^13/ 



-10 meV. The AFM SE energies 
iV[A^T(AcT + Ul')] ^ -2 meV and 
/{U2 + 6J) ~ -0.6 meV. Here iis ^ 0.08 
eV, U^' = [/ -f 2J ~ 6.6 eV and Act ~ 1 eV is the 
charge-transfer energy to put an electron from a Cr site 
to a Ti site, as estimated from Fig . [Ijf). Hopping pa- 
rameters are scaled from Ref. |28[ through the law of 
Ref. (35I I. inversely proportional to the fifth-power inter- 
distance. There are two possible contributions for Epf^ 
(the two opposite paths in the plane of Fig. [3]), one for 
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-^AFM ^^'l fo^' ^AFM ^^'i therefore the SE FM solution 
is stabilized by: ^E^°^ ~ 20-2-1 = 17 meV, that 
overcomes the LSDA-|-?7 AFM solution. We stress again 
the importance of the NNN condition in order to stabi- 
lize ferromagnetism: if the two Cr ions were NN, as in 
rutile, they would have been characterized by a strong 
AFM SE. This calculation shows that, although the F- 
center exchange is not active, Cr'^+ impurities contribute 
to magnetism via SE interactions, in keeping with exper- 
iments 3, 16, 17 1 and contrary to pure DFT results nOl, 
or LSDA-I-U calculations (this work). 

It is important here to highlight that the FM SE sta- 
bilization energy strongly depends on t, ^13 and Act 
parameters. In particular, increasing ^13 and decreas- 
ing Act further stabilizes Epf^. The opposite is true if 
ti3 decreases and Act increase. All this suggests that, 
though they are not needed in our SE mechanism for 
FM, structural defects might change magnetic proper- 
ties by changing the local values of ^13 and Act- In 
Ref. [3, Il7{ it is clearly stated that fast-grown sam- 
ples, structurally defective, show a robust FM, contrary 
to non-defective samples. We suggest therefore to con- 
firm or reject our proposed scenario by polarized EXAFS 
measurements with micro-beam at Cr K-edge, that can 
analyze the local environment of Cr atoms thereby show- 
ing whether ^13, the most critical of the three parameters, 
increases or decreases in fast-grown samples. 




= oxygen vacancy 



FIG. 3: (Color online) Schematic description of SE paths 
around the oxygen vacancy for charge-compensated anatase 
(see text). Ovals represent dzx orbitals. 

Finally, in the absence of O vacancies (Cr'*+), the pre- 
vious mechanism is still active, but in competition with 
another SE interaction, mediated by the bridging O, 
where 2 electrons from the same 2p-orbital at the bridg- 
ing 0-site move toward eg empty orbitals of Cr'*^" ions. 
This interaction is AFM because of the Hund's coupling 
with the underlying t2g orbitals, but its magnitude, for 
realistic values of the parameters, does not exceed 2 meV. 
Overall the LSDA-I-C7 FM solution is again stabilized. 

GoWo@LSDA+U results. We double-check our zeroth- 



order LSDA+U calculations through their GW correc- 
tions (at the Go Wo level) by comparing the quasipar- 
ticle GoWo@LSI)A+U DOS with available photoemis- 
sion experiments. As a general trend. Go Wo corrections 
strongly influence LSDA-f [/ DOS, with different contri- 
butions for spin-up and spin-down states. New empty 
impurity states appear, either isolated [Fig. 1(a), 1(b), 
1(f)] or merged with the bottom of the conduction band 
(BCB) [Fig. 1(c), 1(d)]. 

Our Go Wo results in samples with O vacancies (low- 
est LSDA+U formation energy) provide filled t2g impu- 
rity states at 0.20 eV above the top of the valence band 
(TVB) for rutile and at 0.70 eV for anatase (FM solution, 
stabilized by SE) . This allows us explaining with satisfy- 
ing precision photoemission experiments |14| that found 
a state 0.3 eV (1.0 eV) from the TVB for doped rutile 
(anatase), thereby validating our approach. We remark 
that for anatase, had we evaluated the GqWq corrections 
from the AFM LSDA-t-?/ ground state (i.e., without SE 
corrections), we would have found two impurity states, 
at 1.30 and 2.00 eV, against experimental evidence 14|. 

A direct comparison of calculated band gaps with op- 
tical experiments is in principle invalidated by the ne- 
glect of excitonic effects in our approach. However, it 
is interesting to notice that calculated trends for GqWo 
band gaps in doped and undoped rutile are in very good 
agreement with the experimental trend. In particular, 
we observe theoretically a big decrease in the band-gap 
of doped rutile (2.00-1.90 eV), compared to the undoped 
sample (3.70 eV), in keeping with the experimental trend 
of 1.98 eV for a doped sample [si'], 3.21 eV for the un- 
doped one [13] • 

Conclusions. Our results show that the central role 
in both the induced magnetism and the position of im- 
purity states is covered by the filling of the t2g levels of 
the dopant. In fact, GW@LSDA+[/ DOS shows that 
the band gap is determined by t2g levels of the dopant 
and allows us to exclude the presence of impurity levels 
at the oxygen vacancy. The O vacancy is a F+^-center 
and the associated absence of electrons also implies that 
the F-center exchange, usually invoked to stabilize FM 
in these systems [^, Il4l ] , can not be active in Cr-doped 
Ti02. The role of t2g-electrons for magnetism is high- 
lighted by the new FM stabilization mechanism that we 
propose for anatase, based on the SE path of Fig. [31 
Such a mechanism might be favored by oxygen vacan- 
cies, as the presence of a bridging oxygen opens a com- 
peting AFM SE channel, but it is also extremely sensitive 
to local structural distortions through hopping integrals. 
Actually, the relation of hopping and charge-transfer pa- 
rameters with structural degrees of freedom (vacancies, 
structural defects, bond distance) becomes of crucial im- 
portance in such an approach and can indeed explain the 
enormous variety of results on magnetic ground states 
of a wide class of magnetically doped oxides. We sug- 
gested an experiment that can in principle confirm or re- 
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ject our model. Methodologically, the present approach 
with LSDA+C/ and SE corrections represents a general 
scheme that could be applied to other TMO in order to 
characterize magnetic ground states. In parallel, GqWq 
calculations on top of LSDA+?7(+SE) corrected results 
allow an independent check of the LSDA+C/ ground state 
through a direct comparison of the quasiparticle spec- 
trum with photoemission experiments. 
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